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New ferromagnetic nitrides, CaN and SrN, and their synthesis process 
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We introduce new type of ferromagnets, CaN and SrN, which were designed using first-principles 
calculations. These are half-metallic ferromagnets and they have magnetic moments of 1 fiB per 
chemical formula unit. Out of the typical structures of binary compounds, the rock-salt structure 
is the most stable form for both CaN and SrN. The majority of the magnetic moment of these 
compound originates from the N sites since the p states of N are spin-polarized. Their formation 
energies were calculated and the results show that it should be feasible to synthesize these materials. 
The structural stability of CaN was confirmed by performing first-principles molecular dynamics 
simulations. We propose a synthesis process for CaN basd on the first-principles. 

PACS numbers: 71.20.Dg 71.20.-b 75.50.Cc 75.50.-y 
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Synthesis of substances which do not exist in nature 
promises to have an enormous potential to produce ma- 
terials having excellent properties. Recent experimental 
reports of the syntheses of PtN,[3 IrN, IrNa and OsNzJl 
demonstrate that it is possible to synthesize new materi- 
als using special experimental techniques. 

On the other hand, the reliability of standard first- 
principles calculations has been widely recognized. They 
are very powerful tools for predicting the physical and 
chemical properties of materials including the individ- 
ual characteristics of elements based on their electronic 
configurations. In particular, structural properties and 
physical quantities calculated from them highly accu- 
rate are in good agreement with experimentally mea- 
sured values. To date there has been no example of an 
incorrect prediction of the magnetic properties of sim- 
ple transition metals compounds using a spin-dependent 
generalized gradient approximation (spin-GGA) for an 
exchange-correlation functional. These facts imply that 
it is possible to obtain accurate predictions of struc- 
tural and magnetic properties without performing exper- 
iments. This allows us to design new materials before 
they are synthesized in the laboratory. Our aim is to de- 
sign new materials and discover the synthesis processes 
of these new materials using theoretical techniques. 

In this study, we focus on nitride compounds. It 
is well known that nitrides prefer the ferromagnetic 
state. For example, most rare earth nitrides become 
ferromagnetic 13] and usually crystallize in a rock-salt (RS) 
structure. These facts encourage us to consider replacing 
the rare earth elements with Ca and Sr since their ionic 
radii are close to those of the rare earth ions, and from our 
previous results we confirmed that Ca- and Sr-pnictides 
are ferromagnetic^ . Such replacement is natural in view 
of the Zintl phasej^- There is thought to a probability of 
being able to synthesize RS-CaN and RS-SrN if appropri- 
ate experimental techniques are employed. In fact, PtN 
has been synthesized in a high pressure environment [ij . 

Many predictions have been madejg, S ll> 1^ based 
on density functional theory (DFT)[l^ using the lo- 
cal spin density approximation (LSDA) or the spin- 



GGA, 



and these predictions have been experimentally 



i tnese predictions na 
confirmed. IH El [la, ii , J^ JJ] Spin-GGA calculations 
of the electronic structures and the magnetic properties 
of materials which consist of atoms which have only s 
and p electrons tend to be the most reliable, and those 
of second-row elements such as carbon and oxygen, give 
satisfactory results. 

We performed first-principles calculations based on 
DPT using the LSDA or the spin-GGA for RS-CaN and 
RS-SrN. Typical DPT codes were used to confirm our 
results |l3, [H, lly. The total energies of some struc- 
tures and each magnetic state and detailed electronic 
structures (e.g. densities of states (DOS) and magnetic 
moments) were calculated by the APW+lo method. |17| 
The value of RMrKmax was fixed at 8.00, where Rmt 
is the minimum muffin-tin (MT) radius and K„iax is 
the maximum reciprocal lattice vector. The MT radii 
of Ca and Sr (MTca(Sr))are 0.22a, 0.27a and 0.45a for 
the zinc-blende(ZB), RS, and CsCl structures, respec- 
tively, where a is the lattice constant. The MT radii 
of N (MTn) are 0.18a, 0.221a, and 0.3681a for the ZB, 
RS, and CsCl structures, respectively. Por the wurtzite 
structure MTca(Si) and MTn were set to fixed values 
for each volume and ratio between MT^afSr) and MTn 
is 1.3:1 and the sum of both MT radii is about 95% of 
atomic distances. We used an angular momentum expan- 
sion up to Imax = 10. The energy convergence criterion 
was set to 0.001 mRy. Twenty k-points were taken in the 
irreducible Brillouin zone. To confirm magnetism in the 
four types of structures, we performed more accurate cal- 
culations which employed ten times larger k-point sam- 
pling in the irreducible zone. By doing this we were able 
to confirm the definite existence of magnetism for these 
cases. To check the magnetism, we performed calcula- 
tions for RS-CaN using the KKR-CPA-LDA method Ts]. 
Por Morruzi-Janak- Williams method was used to calcu- 
late the exchange-correlation functional 20]. Por this 
calculation we used an angular momentum expansion up 
to Imax = 2. We used the Espresso-code 19J to perform 
first-principles molecular dynamics simulations and also 
to check the magnetism of RS-CaN . The energy cut-off 
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FIG. 1: (Color online) Total energies and magnetic moments 
of (a) CaN and (b) SrN are plotted as a function of volume. 
The magnetic moment per chemical formula unit is shown on 
the right vertical axis. 
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FIG. 2: (Color online) The DOSs of RS-CaN. The partial 
DOS is defined in the muffin-tin sphere used in the band cal- 
culations. 



of plane waves was 80 Ry. 8x8x8 k-point meshes were 
used for first Brillouin zone integration. Other details 
concerning the calculations are described later. 

In Figs, n^a) and^b), the total energies of CaN and 
SrN are plotted as a function of volume per chemical for- 
mula unit. We chose the RS structure along with the 



TABLE I: 


The magnetic moment of CaN. In units of /ib- 


structure 


RS ZB CsCl 
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1.000 


1.000 


Ca 


0.038 


0.038 


0.033 
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0.794 


0.787 


0.822 


Interstitail 


0.168 


0.175 


0.145 



other typical structures of binary compounds, namely, 
ZB, wurtzite and CsCl structures. The most lowest struc- 
ture is the RS structure, while the wurtzite structure is 
the second lowest stable structure. The wurtzite struc- 
ture was optimized with respect to both c/a and the 
internal parameter. We checked the stabilities against 
tetragonal distortion for the RS, ZB and CsCl struc- 
tures and confirmed that the RS and the CsCl struc- 
tures are stable. In particular, the stability of the cu- 
bic structure for the RS structure is doubly confirmed 
on the basis of these total energy calculations and the 
first-principles molecular dynamics simulations described 
below. By contrast, tetragonal distortion reduces the to- 
tal energy of the ZB structure, which is similar to the 
case of CaAs in the ZB structure^. In the ZB struc- 
ture with tetragonal distortion, c/a=1.0 is the maximum 
energy point and c/a=0.6 is the minimum energy point. 
The energy difference between these two points is about 
0.53eV(0.039Ry) for CaN and 0.23eV (0.017Ry)for SrN. 
Based on these differences, the distorted ZB structure is 
thought to have almost the same energy as the wurtzite 
one (Figs, da), (b)). 

The four structures used in the calculations for CaN 
and SrN are of half- metal as well as ferromagnetic in spite 
of not containing any transition or rare earth metals. In 
other words, CaN and SrN are ferromagnetic without 
having any d or / electrons. The magnetic moments are 
plotted in Figs. Q^a) and (b). The calculated magnetic 
moments for both CaN and SrN are 1 /zb, which is evi- 
dence of their half-metallicity. The DOSs of RS-CaN are 
shown in Fig. |21 The total DOS indicates completely 
spin-polarized electronic states at the Fermi level. In 
the other structures, CaN and SrN are also half-metallic 
as shown in Fig. ^ while the magnetic moments of 
Ca(P,As,Sb) and Sr(P,As,Sb) depend on the structures 
and lattice constants or volumes. |j| The partial DOSs 
indicate that the principal component at the Fermi level 
consists of p states of N. As expected, the majority of the 
magnetic moment resides on the N atom as shown by Ta- 
ble^ The magnetic moment at the Ca or Sr sites is quite 
small being less than that of the interstital region. 

We show in Fig. O the difference in the total ener- 
gies for the ferromagnetic, antiferromagnetic and non- 
magnetic states for RS-CaN. The energy difference be- 
tween the ferromagnetic and antiferromagnetic states for 
RS-CaN is more than 0.1 eV. These values are quite large. 

We further confirmed the structural stability of CaN 
having the RS structure by using first-principles molec- 
ular dynamics simulations. The simulations were per- 
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FIG. 3: (Color online) Total energies as a function of the 
lattice constant are compared for the ferromagnetic (FM), 
antiferromagnetic (AFM) and nonmagnetic (NM) states of 
CaN, where RS structure is assumed for both substances. 
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FIG. 5: (Color onhne) The enthalpy of /3-Ca3N2 and the 
sum, 2CaN+Ca, as a function of pressure. Three types of 
structure of Ca are assumed and plotted. 
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would be stable if it can be synthesized. 

We estimated the formation energy defined by 
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FIG. 4: (Color online) A promising synthesis process for 
CaN. 



formed using the cubic unit cell, which contains eight 
atoms and is four times larger than the primitive cell, 
as the starting cell. We did not assume any symmetry 
for the atomic positions, that is, we employed a triclinic 
structure. We performed two types of molecular dynam- 
ics simulations; (1) a first-principles damped dynamics 
simulation using a variable cell, and (2) a first-principles 
constant-pressure and temperature molecular dynamics 
simulation using the Parrinello-Rahman method, where 
we used the velocity scaling method to control the tem- 
perature. The purpose of simulation (1) was to confirm 
whether the RS structure is at least metastable with re- 
spect to changes in the internal atomic positions and in 
the shape of the unit cell, while the purpose of simula- 
tion (2) was to confirm the stability of the RS structure 
at finite temperature and atmospheric pressure. From 
simulation (1) we confirmed that atomic forces become 
negligibly small (less than 1.0 x 10~^ Ry/a.u.) and 
that atomic positions in the RS structure are at least 
metastable. From simulation (2) we confirmed that the 
RS structure is stable at atmospheric pressure and finite 
temperatures (130 K, 210 K, and 760 K). These results 
support the possibility that CaN having the RS structure 



^tot - (^Ca(Sr) + EnJ2) for the RS, ZB 
'ca(Sr) is the total energy of the 
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and CsCl structures. Et 

bulk Ca (Sr) per atom at a atmospheric pressure which 
we evaluated using the fee structure. E-^^ denotes the 
total energy of an isolated N2 molecule. The energy of 
the molecule is calculated using a large unit cell (a=18 
A, for an fee cell cotaining one molecule). The calculated 
formation energies are -11.2 eV(RS), -10.7 eV(ZB), and 
-10.8 eV(CsCl) for CaN, and -10.8 eV(RS), -10.5 eV(ZB), 
and -10.5 eV(CsCl) for SrN, respectively. These results 
suggest that CaN and SrN having the three structures 
investigated, will not decompose once they have been 
formed. Out of these structures, the formation energy 
of the RS is lower than that of the other structures for 
both CaN and SrN. 

Under normal conditions, Ca and N usually form 
Ca3N2 in the /a3(a-Ca3N2) structure. This structure 
is transformed to /?-Ca3N2 whose structure is P3ml|2l| 
when the temperature is increased. We compared the for- 
mation energies of RS-CaN with that of Ca3N2 for both 
a- and /9-Ca3N2. The formation energy of the a phase 
is about 2.9 eV lower than that of 2 CaN -I- Ca. The 
formation energy of the /? phase, however, is about 2.5 
eV higher than that of 2 CaN + Ca. This fact implies 
that if CaN is formed using an appropriate method, it 
may be more stable than /3-Ca3N2. 

We now proceed to propose ideas for how CaN and SrN 
can be synthesized in the laboratory. In this paper, we 
only consider RS-CaN. Based on the above results, sev- 
eral possible methods and processes suggest themselves. 
One of them is high pressure and high temperature syn- 
thesis. We propose the synthesis process based on first- 
principles calculations. A scenario for synthesizing RS- 
CaN is illustrated in Fig. 0] Firstly, heat a-Ca3N2 until it 
is transformed into /3-Ca3N2 whose structure is P3to1|21J 
and then compress it using a pressure higher than 50 GPa 
until the reaction /?-Ca3N2 -^ 2CaN -I- Ca occurs. After 



that, cool the product and decompress it to atmospheric 
pressure. In order to evaluate the transition pressure, we 
give the enthalpy differences in Fig. EJ where we show 
the enthalpy of /?-Ca3N2 and that of the sum of the two 
phases, 2CaN+Ca. We assumed the structures of Ca as 
fee, bcc, or sc. As the pressure is increased, the structure 
of Ca is transformed from the fee into bcc (at 20 GPa) 
into sc (at 32 GPa). The enthalpy of 2CaN+Ca (sc) 
becomes lower than that of (3-03.3^2 at about 50 GPa. 
Based on the calculations of the formation energy of RS- 
CaN and the constant pressure and temperature first- 
principles molecular dynamics simulation, 2CaN+Ca(sc) 
should not decompose once it has been formed. This re- 
sult strongly suggests that /3-Ca3N2 is transformed into 
RS-CaN with sc Ca. The transition temperature from 
a-Ca3N2 to the /3-Ca3N2 is only represented by a dashed 
line, since it has not been precisely determined experi- 
mentally. Based on this scheme, laser heating may heat 
the sample to a temperature that is too high for the /?- 
phase. Accurate experimental determination of the phase 
diagram of Ca3N2 is required to determine the transition 
temperatures. 

In the view of experiment, we might consider /3-Ca3N2 
-f N2 -^ 6 CaN reaction. In this process N2 is liquid and is 
used to compress hydrostatistically. We have not checked 
the influence of pressure for that process yet, and we do 
not deny the possibility. 

The other method having a high probability of suc- 
cess is molecular beam epitaxy. To synthesize the RS 
structure, we propose a synthesis process using a su- 
perlattice structure, similar to the fabrication process of 
CrAs/GaAs which used a superlattice structure "111. On 
surfaces, processes may occur which are very different 
from those that occur in the bulk material and, conse- 



quently structures that are different from those in the 
bulk may be produced. This is the case of ZB-CrAs, 
which is unstable according to total energy calculations. 
We do not deny the existence of other synthesis processes, 
for example, heating to a liquid state and then compress- 
ing to an appropriate pressure, similar to the process for 
synthesizing PtN. Our suggestions are only two of several 
possible processes. We hope that our suggestion inspires 
experimentalists to synthesize CaN and SrN. 

In conclusion, we have designed and investigated new 
ferromagnetic nitrides, CaN and SrN. These are half- 
metallic ferromagnets as a result of the polarization of 
p-orbitals of N. The process for synthesizing RS-CaN us- 
ing a high pressure is proposed based on the results of 
the first-principles calculations. The development of ex- 
perimental techniques may allow us to synthesize such 
new materials. 
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